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Precision computation of the leptonic D s -meson decay constant in quenched 
QCD 
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Department of Physics, Humboldt- University Berlin 

We summarize a computation of the leptonic decay constant fb a of the D s -meson in quenched QCD on the lattice. We 
perform a direct simulation at the masses of the strange and the charm quarks at four different lattice spacings from 
approximately 0.1 fm to 0.05 fm. Fully non perturbative 0(a)-improvement is employed. After taking the continuum limit 
we arrive at a value of fb s = 252(9) MeV, when setting the scale with the Kaon decay constant Fk ~ 160 MeV. Setting 
the scale with the nucleon mass instead leads to a decrease of about 20 MeV of F&. . 



1 Introduction 

To get reliable estimates of weak decay constants like 
Fb lattice QCD has often to be supplemented by chi- 
ral extrapolations and/or heavy quark effective theory. 
These introduce substantial systematic errors In 
addition, the usual lattice errors like statistical errors, 
discretization errors, finite volume effects, contamina- 
tion from excited states, perturbative renormalization 
and quenching [ |2] have to be understood or elimi- 
nated. 

The D s -meson however is special in this context. It 
consists of a strange and a charm quark. Both can 
be implemented precisely directly on the lattice as has 
been done recently by the ALPHA collaboration [|2||2]. 
Thus neither a chiral extrapolation nor an extrapola- 
tion to heavy quarks (or similar strategies) have to be 
used. Still we expect the D s -meson to be similar to 
other heavy light systems. Therefore it can be studied 
to understand all the other error sources mentioned 
above. 

The goal of our work, which has been recently pub- 
lished in [ E| , is a precision computation of i*b B in 
quenched QCD. We aim at a combined error of three 
percent. To eliminate the discretization errors, we 
perform fully non perturbative 0(a)-improvemcnt and 
simulate at four different lattice spacings. This allows 
us to take a reliable continuum limit. Finite volume 
effects have been shown to be negligible in [ . At 
our masses they are expected to be even smaller. We 
define our plateaus such that the contamination by ex- 
cited states stays below five per mille. Finally all the 
uncertainty due to perturbative renormalization has 
been eliminated by using the non perturbative renor- 
malization techniques of the ALPHA collaboration [ 



El ■ The only systematic uncertainty we cannot deal 
with at present is the quenching error. 

A precise value of i<b a in quenched QCD is desirable 
since together with a computation of the decay con- 
stant in the static approximation it will enable us to 
see how far the heavy quark effective theory can be ap- 
plied safely. In unquenched simulations we can then 
rely on this experience in the computation of Fg ■ 

The D s -meson is stable in QCD. It decays weakly by 
an emission of a W-boson into a lepton and a neutrino. 
The branching ratios can be measured experimentally. 
They are summarized in Given the CKM-matrix 
element V^ s , Frj B can then be determined, since the 
standard model predicts 



BR(D S -> bv) 

fV 

8tt s 5 



(1) 
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Recent experimental data are shown in table G] [ E] • 
The status of lattice computations was reviewed in [ 



F Ds [MeV] 



ALEPH 
DELPHI 
L3 
CLEO 
BEATRICE 
E653 



285 ± 19 ± 40 

330 ± 95 
309 ± 58 ± 50 
280 ± 17 ±42 
323 ± 44 ± 36 
194 ± 35 ± 24 



"This work is part of the ALPHA research programme. JR 
thanks the organizers of the conference for the kind invitation. 
tHU-EP-03/26, SFB/CPP-03-09 



Table 1. Experimental data for Fb B 



10 with a quenched world average of Fq b = 230 ± 
14 MeV. Table summarizes some of the latest re- 
sults. 
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Fp 3 [MeV] 

UKQCD [HU 229(3) + 23- 12 
APE [El 234(9) + 5-0 
MILC [HS| 223(5) + 19 - 17 

Table 2. Lattice results for i<b B in quenched QCD. 

2 Strategy 

The decay constant i<b B is defined by the QCD matrix 
element 

(0|A M (0)|D s (p))=ip M F Ds (2) 

of the axial current = sj^j^c. To formulate this 
problem on the lattice wc eliminate the bare param- 
eters of the QCD Lagragian in favour of physical ob- 
servables in one chosen hadronic scheme. Our strat- 
egy [ El E3 E] is to use the kaon decay constant F K 
to set the scale, that means to compute the lattice 
spacing a in physical units as a function of the bare 
coupling g . The bare strange quark mass and the 
bare charm quark mass are eliminated by the masses 
rriK and mD, of the kaon and the D s -meson, respec- 
tively. We neglect isospin breaking and take the quark 
mass ratio M s /Mii g ht = 24.4 ± 1.5 from chiral pertur- 
bation theory [ITo]. Here Mi ight = ±(M U + M d ). The 
capital letters denote the renormalization group in- 
variant quark masses, which are somewhat larger than 
the MS-masses. In quenched QCD they take the val- 
ues M s = 138(5) MeV [Ej and M c = 1.654(45) GeV [ 
0] . The results mentioned above have been obtained in 
units of the Sommer scale tq which is derived from the 
force between static color sources, r^/a has been com- 
puted to a high precision in [ 1161 117| as a function of the 
bare couling. The Sommer scale is affected by lattice 
artifacts only at 0(a 2 ) in the quenched theory. Setting 
the scale with Fk — 160 MeV is equivalent to using 
ro = 0.5 fm while the nucleon mass tojv = 938 MeV 
roughly corresponds to ro = 0.55 fm [EJ- This devia- 
tion is the typical size of the quenched scale amiguity. 

3 Theory 

We perform numerical simulations in 0(a)-improved 
lattice QCD using Schrodinger functional boundary 
conditions [ on a L 3 x T space-time cylinder. 

For unexplained notation we refer to [EOj- We define 
the meson sources 

= a 6 J2Uy)^Q^), (3) 

°' = a6 E^(y)75caz), (4) 

with flavour indices i and j at the Xq = and x = T 
boundary, respectively. From these we compute the 



correlation functions 

/a(*o) = -\{OA Q {x)), (5) 

h = -^(°'°)- ( 6 ) 

Here A^(x) denotes the improved axial current. It 
receives a scale independent multiplicative renormal- 
ization Z A on the lattice. In terms of these correlation 
functions -Fb B can be written as 

Fd, = -2Z A (l + b A (am qyi + am q j)/2) 

x /aM ( TODB L 3 )- 1 / 2 e^- T / 2 ) mD = 
v/i 

x{l-^e^ A -77°e-( T ^> m «} 
+0(a 2 ). (7) 

Here the factor (mD s L 3 ) -1 / 2 takes into account the 
normalization of one particle states. The contribu- 
tion f-y cancels out the dependence on the meson 
sources. Because of the exponential decay of the cor- 
relation function /j[ the product in Q is expected 
to exhibit a plateau at intermediate times when the 
contribution r] ! ^ s e^ XoA of the first excited state and 
the contribution rp A e~^ T ~ x °^ mG from the ++ glueball 
both are small. A plateau average can then be per- 
formed to increase the signal and is understood in Ql . 
Further explanations for equation J7J and details can 
be found in [ I21j . 

4 Parameters 

We discretize the space-time cylinder using four differ- 
ent lattice spacings a but keep L and T = 2L approx- 
imately constant in physical units. To this end we use 
the same bare couplings that have been used in the 
determinations of the strange and the charm quark 
masses in [ El E] by the ALPHA collaboration. From 
this work we also take the hopping parameters for the 
quarks. Our choice of parameters is shown in table El 
In table 0] we show that with this choice indeed the 






^mcas 


L/a 


K s 




6.0 


380 


16 


0.133929 


0.119053 


6.1 


301 


24 


0.134439 


0.122490 


6.2 


251 


24 


0.134832 


0.124637 


6.45 


289 


32 


0.135124 


0.128131 



Table 3. Statistics and parameters for our simulations. 



physical conditions are constant to a sufficient preci- 
sion. 
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a 


L/r Q 


romD B 


6.0 


2.98 


4.972(22) 


6.1 


3.79 


4.981(23) 


6.2 


3.26 


5.000(25) 


6.45 


3.06 


5.042(29) 



Table 4. Demonstration of constant physical conditions. 








6.0 


0.540(14) 


6.1 


0.576(13) 


6.2 


0.598(16) 


6.45 


0.614(15) 


c.l. 


0.638(24) 



Table 5. Simulations results and continuum limit for Fd s - 
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Figure 1. Behaviour of aF^ TC . The full symbols denote 
the plateau range. 



Continuum extrapolation of r F D 
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Figure 2. Continuum extrapolation of fb s . 

Alamos group [ I23|. Since this involves an extrapo- 
lation of their data we have also used 1-loop pertur- 
bation theory Then we get r F Bs = 0.631(24), 
which is in perfect agreement with our main result, 
r F Ds = 0.638(24), or, using r = 0.5 fm, F Ds = 
252(9) MeV. The continuum extrapolation is shown 
in figure |21 



5 Computation of the decay constant 

To compute the decay constant Fu s we use the com- 
bination of correlation functions. For all parame- 
ter choices we find plateaus as functions of xq. These 
plateaus are shown in figure^ Their extent is roughly 
from 4ro to 5ro. At small respectively large times we 
fit to the functions Fo s (xq) the expected contributions 
of the first excited state and the glueball Q . For de- 
tails see [ [S] . We define the plateau such that their 
sum stays below 5 per mille. 

Since we deal with heavy quark propagators on 
APE 1000 in single precision we have to check that the 
rounding errors are small enough. Our check against 
runs with a double precision code reveals that the 
rounding errors are smaller than one per mille. 

After averaging Fd s over the plateaus defined above 
we get the values shown in table |5] These data can 
be extrapolated to the continuum limit. Here we leave 
out the coarsest lattice. The extrapolation can be per- 
formed in (a/r ) 2 since we employ non perturbative 
O(a) improvement. Here we take 6a from the Los 



6 The quenched scale ambiguity 

To estimate the quenched scale ambiguity of i*b B un- 
der a scale shift of 10 percent, which is typical for the 
quenched approximation, we consider roi*b s = f{z) 
as a function of the meson mass z = r^m^^. We 
expand f(z) around the physical value zq — 4.988 
up to first order. A 10 percent increase of ro cor- 
responds to z — zq = 0.5. With an estimate of 
f'(zo) from a linear fit of roFv a around mo s we get 
f(z) — f(zo) w 0.008(3). Converting back to physi- 
cal units, now using ro = 0.55 fm we find that Fd s 
decreases by 20 MeV, corresponding to eight percent. 
The estimate of f'(zo) is possible since in addition to 
the hopping parameters already discussed we have also 
performed simulations around the charm quark mass 
value. 

7 Conclusion 

The leptonic D s -meson decays can be studied on the 
lattice without chiral extrapolations or heavy quark 
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effective theory. This has enabled us to perform a 
computation of Fd s with a precision that matches the 
precision goals at future experiments, for example, 
CLEO. The precise value of Fd b in quenched QCD, 
together with new precise data in the static approxi- 
mation, will show how far heavy quark effective theory 
can be applied safely. This is of importance for the un- 
quenched computation of in the future. 

We will supplement this analysis with more data 
around the charm mass. Part of this data has already 
been used to estimate the quenched scale ambiguity of 
Fd b under a scale shift of 10 percent. 
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